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ABSTRACT 

Aims. We examine the dependence of the wind-wind collision and subsequent X-ray emission from the massive WR+O star binary 
WR 22 on the acceleration of the stellar winds, radiative cooling, and orbital motion. 

Methods. Three dimensional (3D) adaptive-mesh refinement (AMR) simulations are presented that include radiative driving, gravity, 
optically-thin radiative cooling, and orbital motion. Simulations were performed with instantaneously accelerated and radiatively 
driven stellar winds. Radiative transfer calculations were performed on the simulation output to generate synthetic X-ray data, which 
are used to conduct a detailed comparison against observations. 

Results. When instantaneously accelerated stellar winds are adopted in the simulation, a stable wind-wind collision region (WCR) 
is established at all orbital phases. In contrast, when the stellar winds are radiatively driven, and thus the acceleration regions of 
the winds are accounted for, the WCR is far more unstable. As the stars approach periastron, the ram pressure of the WR's wind 
overwhelms the O star's and, following a significant disruption of the shocks by non-linear thin- shell instabilities (NTSIs), the WCR 
collapses onto the O star. X-ray calculations reveal that when a stable WCR exists the models over-predict the observed X-ray flux by 
more than two orders of magnitude. The collapse of the WCR onto the O star substantially reduces the discrepancy in the 2 - 10 keV 
flux to a factor of 6 at = 0.994. However, the observed spectrum is not well matched by the models. 

Conclusions. We conclude that the agreement between the models and observations could be improved by increasing the ratio of the 
mass-loss rates in favour of the WR star to the extent that a normal wind ram pressure balance does not occur at any orbital phase, 
potentially leading to a sustained collapse of the WCR onto the O star. Radiative braking may then play a significant role for the WCR 
dynamics and resulting X-ray emission. 

Key words. Stars: winds, outflows - Stars: early-type - Stars: individual (WR 22) - Stars: massive - X-rays: binaries - 
Hydrodynamics 

1, Introduction The X-ray emission from a colliding winds binary sys- 

tem acts as a direct observational probe of the postshock 
Hot luminous massive stars possess powe rful radiatively dr iven winds, and hence an indirect probe of the preshock winds (e.g. 
stellar winds (for a recent review see |Puls et al.| EQQSb- In [Stevens et al] Il996l: iPittard & Corcoranl l2002l IpTBecker et al.l 
a binary system consisting of two such stars, the collision |2006l) . R ecently, an analysis o f XMM- Newto n observations of 
of the winds generates a region of high temperature (T > WR22 bv lGosset et al. (20091 (hereafter lG09h characterised the 
10 7 K) plasma which emits prolificalry at X-ray wavelengths X-ray emission using a two-component spectrum consisting of 
(Prilutskii & Usov[ 1976t [Cherepashchug|l976|; [Euo et al.||1990fc a soft component at ~ 0.6 keV and a harder component at 



IStevens et al.lll992l). Depending on the parameters of the winds -2-4.5 keV. However, difficulties were encountered as wind- 

and the orbit, the dynami cs of the postshock gas in the WCR wul d collision models were found to over-predict the observed 

can cover a diverse range (|Stevens et aL||l992|). For instance, in X-ray flux by more than two orders of magnitude, 
long-period binaries (i.e. of the order of years) the postshock 

gas is expected to be quasi-adiabatic for the most part, whereas Considering the separation of the stars in WR22 (Table 

in short-period (i.e. a few days) systems the postshock gas is ex- and the dominant WR wind ram pressure characteristic 

pected to be highly radiative. As a result, eccentric intermediate oj_WR+0 binary systems (e^ WR140 |Williams et alJ 

period systems (of the order of 100's of days), such as WR22 LL222t |Zhekov & Skinned |2000t Pollock et al.| |2005t 

(e * 0.559, P * 80 days - see Tables H] and 0), provide the IPittard & Dougherty! |2006|), the WCR will reside in the 

prospect of transitioning between these two extremes. Such a wind acceleration region of the O star's wind throughout the 

transition in the state of the postshock gas has also been found orbit - Consequently, lower preshock velocities will increase 

in simulation s of shorter period eccentric OB star systems by the importance of radiative cooling in the postshock gas, 

iPittardl Jl998h and lPittardl d2009b . affecting the stability of the WCR. The interplay between 

the stellar radiation fields may also significantly reduce the 

acquired preshock velo cities through r adiative inh i bition 

* F. R. S.-FNRS (Belgium) (IStevens & Pollockl 1 19941) and/or braking dGavlev et al.l 1 19971 
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Table 1. Adopt ed system pa r amete rs for WR 2 2. R eferences are 
as foll ows: 1 = iRauw et alJ (Il996h 2 = |G09, 3 = iBohlin et al.1 
(119781) . 4 = lDiplas&Savagd d 1994b . 



Parameter 


Value 


Reference 


Orbital period (d) 


80.325 


1 


a (au) 


1.68 


1 


Eccentricity (e) 


0.559 


1 


Distance (kpc) 


2.7 


2 


ISM column (10 21 cm- 2 ) 


2.5 


3+4 



Table 2. Stellar parameters used to calc ulate the line driving of 
the stellar winds, k and a are t he IC AKl line driving param eters. 
References are as follows: 1 = [van der Huch t et al.l (1 19811). 2 = 
ISchweickhardt et al.l (fl999l). 3 = IRauw et al.l (119961) 4 = iRauwl 



(119971). 5 =lHamann et alJ l|200i 
IVink et alJ d2~ 



. 6 = lMartins et all (120051) . 7 = 
20011) . 8 = ICrowther etaD (ll995bl) . 



WR 



O star 



Parameter 


Value 


Reference 


Value 


Reference 


Spectral Type 


WN7 


1 


09V 


2 


M(M Q ) 


72 


3 


25.7 


3 


A (Ro) 


20 


4 


11 


4 


r eff (K) 


44700 


5 


33000 


6 


L* (1O 6 L ) 


1.44 




0.13 






0.42 




0.30 




Of 


0.47 




0.52 




M(M Q yr 1 ) 


1.6 x 10~ 5 


4 


2.8 x 10~ 7 


7 


Voo (km s -1 ) 


1785 


8 


2100 


7 



see also lOwocki & Gayleyll 19951) . In fact, a stable ram pressure 
balance may not be established between the winds. WR22 - 
one of the most massive Wolf-Rayet stars currently known 
(iRauw et al.ll 19961) - may play host to these interesting phenom- 
ena which have the potential to significantly affect the observed 
X-ray emission. 

The influence of orbital motion on the circumbinary 
medium may also affect the observed emission. The shape of 
the wind-wind collision region (WCR) between the stars is 
largely dependent on the ram pressure of the stellar winds. 
However, orbital m otion introduces curvature to the WCR 
away from the stars dWalder et al.lll999[ iFolini & Walderll2 000 ; 
Wal der & Folinil 120031: iLemaster et al.l 120071: iParkin & Pi ttard 
2008; Okazakietal. 2008; Parkin et al. 2009a; Pittard 2009; 



Par kin et al.l 1201 ll Ivan Marie et al.l 1201 lb . This results in the 
spiral-like stru cture seen in the so-called "pin w heel" nebula (e.g. 
iMonnier et aLlll999b iTuthill et al.ll2006l 120081) . As such, the re- 
sulting highly asymmetric gas distribution introduces a viewing 
angle dependence to the emergent X-ray spectrum. 

In this paper the wind-wind collision in WR22 is investi- 
gated using three dimensional adaptive-mesh refinement (AMR) 
simulations aimed at establishing the importance of the wind ac- 
celeration regions, radiative cooling, and orbital motion on the 
dynamics and resulting X-ray emission. The remainder of this 
paper is structured as follows: the hydrodynamical and X-ray 
calculations are described in § [2j The simulation dynamics, re- 
sulting X-ray emission, and suggested revisions to model param- 
eters are given in § [3j A discussion is given in § HI followed by 
conclusions in § [5] 



2. The model 

2.1. Hydrodynamic modelling 

To model the wind-wind collision we numerically solve the 
time-dependent equations of Eulerian hydrodynamics in 3D 
Cartesian coordinates. The relevant equations for mass, momen- 
tum, and energy conservation are: 



d P „ 
- + V.pv 



0, 



+ V • pvv + VP = pf , 

dt 



dpE 



+ V • [(pE + P)v] 



_P_ 



A(T) + pf • v. 



(1) 

(2) 

(3) 



Here E = e + ||v| 2 , is the total gas energy, e is the specific in- 
ternal energy, v is the gas velocity, p is the mass density, P is the 
pressure, T is the temperature, and mu is the mass of hydrogen. 
We use the ideal gas equation of state, P = (y - l)p6, where the 
adiabatic index y = 5/3. 

The radiative cooling ter m, A(T), is calcu lated from the 
APEC thermal plasma code dSmith et al.l l200ll) distributed in 
XSPEC (v 12.5.1). The temperature of the unshocked winds is 
assumed to be maintained at « 10 4 K via photoionization heat- 
ing by_Jhe_stars : _^FOTjhe WR star w e adopt th e abundances 
of ICrowther et al.l (Il995al) and lG09l (see also lHamann et all 
1 19911). and for the Q s tar solar abundances are assumed 
( Anders & Grevessel [l989). The cooling curves are shown in 
Fig. HI 

The body force per unit mass f acting on each hydrody- 
namic cell is the vector summation of gravitational forces from 
each star, and continuum and line driving forces from the stel- 
lar radiation fields. The calc ulati on of the lin e force has been 
described by iPittardl (120091) and IParkin et al.l (1201 ll) . and we 
refer the reader to these works for further d etails (see also - 
ICranmer & O wockilll995t iGavlev et al.lll997l). In brief, the nu- 
merical sc heme i ncorporates the ICastor. Abbott. & Klein (1975)1 
(hereafter ICAKl) formalism for line driving by evaluati ng the lo- 
cal S obolev optical depth r = <x e v t hp[n • V(n • v)]" 1 (So bolevl 
Il960l) and then calculating the vector radiative force per unit 
mass 



grad = 



(T\- a k 



1(h) 



n V(n v) 



pv t h 



(4) 



where a and k are the standard Ic AKl parameters. <x e is the spe- 
cific electron opacity due to Thomson scattering, and v t h is a 
fiducial thermal velocity calculated for hydrogen. A gaussian in- 
tegration is perfor med to correct the line force for the finite size 
of the stellar disk (ICastodll974l:lPauldrach et al.lll986l) . The line 
driving is set to zero in cells with temperatures above 10 6 K, 
since this plasma is mostly ionized. 

We note that the ICAKl formalism was originally developed 
for O star winds, and as such, certain aspects of WR wind ac- 
celeration are not captured. For instance, the presence of two 
acceleration regions, one in the optically thick inner wind and a 
second in the outer wind, and a shallower velocity profile than 
predi cted by ICAKl (|ffillier & Millet] 1 19991: iGrafener & Hamannl 
120051) . Consequently, the IC AKl formalism will underestimate the 
column densities for lines of sight passing through the inner WR 
wind (R < 100 R*). However, the wind- wind collision generally 
occurs close to the O star and, therefore, the details of the WR 
wind acceleration should not significantly affect the WCR dy- 
namics. Importantly, the ICAKl formalism provides a means for 
incorporating the wind acceleration regions. 
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Fig. 1. Cooling curves calculated for WN7 and solar abundances 
(§P. 



ing a second-derivate truncation error check on the intrinsic 0.5- 
10 keV X-ray flux and sufficiently prominent pixels (£ cr i t = 0.6, 
where £ cr i t is the critical trunc ation error a bove which pixels are 
marked for refinement - see Parkin||201i|) are refined. The pro- 
cess of ray-tracing and refining pixels is repeated until features 
of interest in the image have been captured to an effective resolu- 
tion equivalent to that of the simulation (i.e. 4096 x 4096 pixels). 
For d etails of the AIR technique we refer the reader to Parkin 
d20H . 

To calculate the X-ray emission from the simulation we use 
emissivities for optically thin gas in collisional ionization equi- 
librium obtai ned from look-up t ables calculated from the APEC 
plasma code dSmith et al.l l2001) containing 200 logarithmically 
spaced energy bins in the range 0.1-10 keV, and 101 logarithmi- 
cally spaced temperature bins in the range 10 4 - 10 9 K. When 
calculating the emergent flux we use energ y dependent o pacities 
calculated with version c08.00 of Cloudy (Ferland 2000, see also 
iFerland et al.ll 19981) . The advected scalar is used to separate the 
X-ray emission contributions from each wind. 



2.2. The hydrodynamic code 

The hydrodyna mic equations are solved using v3 .1.1 of the 
FLASH code dFrvxell et alJl2000b iDubev et al. 20 Qgl). This code 
uses the piecewise-parabolic method of Colella & Woodward 
(Il984l) to solve the hydrodynamic equations and operates with 
a block-structured AMR grid (e.g. iBerger & 01iger|ll984l) us- 
ing the PARAMESH package jMacNeice et alJl2000h under the 
message-passing interface (MPI) architecture. The simulation 
domain extends from x = y = +1.2 x 10 14 cm and z = 
(0 - 1.2 x 10 14 ) cm - a symmetry about the xy-plane is used. The 
grid is initialized with xxyxz= 16x16x8 cubic blocks each 
containing 8 3 cells. We allow for 5 nested levels of refinement, 
which results in an effective resolution on the finest grid level of 
x x y x z = 4096 x 4096 x 2048 cells. The refin ement of the 
grid depends on a second-derivative error check dFrvxell et al.l 
2000) on p and the requirement of an effective number of cells 
between the stars to accurately describe the WCR dynamics 
(Parkin et al. 2011). Customized units have been implemented 
into the FLASH code for radiative driving, gravity, orbital mo- 
tion, and radiative cooling f or optically-thin plasma (u sing the 
explicit method described in Strickland & Blond mTl995l) . We in- 
clude an advected scalar variable to track the stellar winds. 

The stellar winds are initiated in the instantaneously accel- 
erated and radiatively driven stellar winds simu lations in two 
slightly di fferent ways and we refer the reader to iPittardl (2009 ) 
and iParkin et al.l (1201 ll) for further details. In short, the instan- 
taneously accelerated winds are initiated into a radius which 
adapts to the simulation resolution (which is dependent on the 
separation of the stars), whereas the radiatively-driven winds are 
initiated into a radius of ~ 1.15 R* around the stars at all orbital 
phases. 

The orbital motion of the stars is calculated in the centre of 
mass frame. When the stars are at apastron the WR and O stars 
are situated on the positive and negative x-axis, respectively. 
The motion of the stars proceeds in an anti-clockwise direction. 

2.3. X-ray emission 

Synthetic X-ray spectra and lightcurves are produced by solving 
the equation of radiative transfer through the simulation domain 
using adaptive image ray-tracing (AIR). An initially low resolu- 
tion image equivalent to that of the base hydrodynamic grid (i.e. 
128 x 128 pixels) is constructed. The image is then scanned us- 



3. Results 

For the purposes of our investigation we have performed two 
simulations: one in which the stellar winds are assumed to be 
instantaneously accelerated at the surface of the star (model A), 
and another in which the winds are radiatively driven (model B). 
Our adopted orbital and stellar parameters are noted in Tables [T] 
and [2] To form a basis for the analysis of the hydrodynamical 
simulations, we first estimate the influence of the wind accel- 
eration regions on the WCR dynamics. The dynamics and X- 
ray emission calculations from models A and B are then pre- 
sented, following which we consider some potential revisions to 
the adopted parameters for WR22. 

3.1. Estimating the impact of wind acceleration 

In model A the stellar winds collide at their terminal velocity 
throughout the orbit, albeit with a small additional velocity com- 
ponent due to contraction/expansion of the binary separation. 
The wind momentum ratio r] - MoVo/^wrVwr = 0.02, and 
the relative distance of the WCR from the stars (along the line 
of centres), stays roughly constanfl Consequently, a stable ram 
pressure balance exists between the winds at all orbital phases. 
The WCR dynamics are la rgely dictated by the importance of 
radia tive cooling (see e.g. IStevens et al.lll992t IParkin & Pittard 
I2010h whic h can be quantified using the cooling parameter 
y (IStevenset al.lll992ln . The postshock WR's wind will cool 
rapidly at all orbital phases (^wr - 0.35 - 1.3 - Fig.0, whereas 
the postshock O star wind will be quasi-adiabatic throughout the 

1 The thermal pressure of postshock gas offsets the position of the 
stagnation point from that attained by the ram pressure of the winds 
alone (e.g. Kenny & Taylor 2005; Gayley 2009). Due to the variation in 
radiative cooling for the postshock gas through the orbit the offset due 
to postshock thermal pressure also varies. 

2 The cooling parameter^ = f C ooi/*flow = Vgdi 2 /M_ 7 , where t coo \ is 
the cooling time, £ flow is the flow time, v 8 is the preshock gas veloc- 
ity (in 10 8 cm s _1 ), d\ 2 is a characteristic distance (in 10 12 cm) taken 
here to be the stagnation point radius, and M_ 7 is the mass-loss rate (in 
10~ 7 M© yr _1 ). x is the ratio of the cooling time to the flow time; x > 1 
indicates quasi-adiabatic gas, whereas x ~ 1 indicates that postshock 
gas cools rapidly. Note that in our calculation of ^ W r we have accounted 
for the emissivity of WN7 abundance gas at 10 7 K being roughly a fac- 
tor of 1.7 higher than that of solar abundance gas (see Fig.[T), hence our 
values of ^ W r are lower than those quoted bv lG09l 
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0.4 0.6 
Orbital Phase 




0.4 0.6 
Orbital Phase 



Fig. 2. Cooling parameter, ^ , calculated for the WR's wind (left panel) and O star's wind (right panel), respectively. Note the 
difference in scale between the plots. 



orbit (xo - 10 - 30). The velocity and density shear across the 
contact discontinuity (CD) will cause Kelvin-Helmholtz (KH) 
instabilities to develop. As the postshock WR wind cools to 
form a thin dense layer, ripples introduced by KH and Rayleigh- 
Taylor ( RT) instabilitie s will seed the growth of thin- shell insta- 
bilities dVishniadl 19831) . However, the hot, quasi-adiabatic post- 
shock O star wind will act as a cushion and prevent the evolu- 
tion of t he bending mod es to non-linear thin- shell instabilities 
(NTSIs - IVishniadll994l) . 

The inclusion of the wind acceleration regions introduces a 
number of important differences to the WCR dynamics. Firstly, 
the separation of the stars is sufficiently small that the winds do 
not reach their terminal velocities before they collide. Therefore, 
the effective wind momentum ratio is a function of orbital phase 
as it depends on the ram pressure balance between the accelerat- 
ing winds. Because the ram pressure of the WR wind dominates, 
the WCR (if one exists) will occur closer to the O star than in the 
instantaneously accelerated winds case. Hence, at phases close 
to periastron the WCR will reside deep in the acceleration region 
of the O star's wind, and the acquired preshock velocities will be 
well below the terminal wind velocity. The impact of wind ac- 
celeration on the importance of cooling in the postshock gas can 
be estimated if one approximates the winds as ft- velocity laws 
(i.e. v(r) = VooQ - RJrf with J3 = if] and solves for the ram 
pressure balance between the winds to determine effective values 
for x- In this case > the postshock O star's wind remains quasi- 
adiabatic for the majority of the orbit (0.08 < (p < 0.92) but 
cooling becomes important at phases close to periastron (Fig.0. 
Furthermore, at cp ^ 0.95 a ram pressure balance is not attained 
between the winds and the WR wind now collides against the 
O star - a corresponding step in ^wr occur s indicating the sud- 
den increase in the stagnation point radius. iGayley et all (1 19971) 
note that for WR22 a WCR may be established due to the ra- 
diative braking of the WR wind by the O star. Evaluating the 
equations of IGayley et al.l (Il997l) . we find that radiative brak- 
ing should occur at periastron and apastron (Fig. [3]- although 
at apastron the WR wind has reached the WCR before the on- 
set of braking). For comparison, hydrodynamic simulations with 



3 Observational and th eoretical studies suggest that J3 > 1 is appro- 
priate for WR wind s (e.g. iKoe nigsberger 1990; Auer & Koenigsberger 
1994; Springmann 1994: iGavlev et all 119951 : iHillier & Milled 11999: 
Grafener & Hamann 2005). Howeve r, as th e numerical scheme for the 
radiatively-driven winds is based on ICAKl adopting J3 = 1 allows pre- 
dictions to be made for the results of model B. 




HS - periastron 
HS - apastron 
G97 - periastron 
G97 - apastron 



200 300 400 

Distance from WR star (Rsol) 



500 



Fig. 3. The WR's wind velocity along the line of centres cal- 
culated using hydrodynamic simulations with radiatively driven 
stellar winds ( HS), and using the analytical radiative braking 
prescription of IGayley et al.l (1 19971) (G97). The crosses indicate 
the position of the O star at periastron (~ 158 R©) and apastron 
(~ 560 R Q ), and the cross width corresponds to the diameter 
of the O star. Note that orbital motion is not included in these 
calculations. 



radiatively driven winds, and separations corresponding to apas- 
tron and periastron, have been performecQ (using the code de- 
scribed in § 0. Although the WR's wind is not as effectively 
halted before it reaches the O star in the hydrodynamic simula- 
tion at periastron separation, the calculations agree in that some 
braking should occur. Note that in these calculations we have 
used the O star's ICAKl parameters (ko = 0.30 and a<j = 0. 52) to 
assess the radiative braking of the WR's wind. If the WR 's ICAKl 
parameters were used (i.e. ko = 0-42 and ao = 0.47) radiative 
braking would be more effective. 

Notwithstanding the possibility of radiative braking, the im- 
portance of cooling for both winds will permit the growth of 
NTSIs close to the WCR apex. The subsequent disruption, and 
possible collapse, of the WCR will significantly affect the dy- 
namics and resulting X-ray emission. 



4 In these simulations the O star radiation field is included but its 
wind is not, hence no wind- wind collision shocks occur. 
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3.2. Model A Dynamics 

The stronger WR wind causes the WCR to form a bow- shock 
which is concave from the perspective of the O star. At apas- 
tron one sees that the influence of orbital motion on the arms 
of the WCR is relatively minor due to the low velocities of the 
stars compared to the wind velocities (Fig. [4]). The postshock O 
star wind is quasi-adiabatic, whereas the postshock WR wind is 
reasonably hot close to the apex of the WCR, but cools rapidly 
downstream to a temperature of ~ 10 4 K. A comparison of the 
density and temperature snapshots in Fig.|4]highlights the region 
of the WCR composed of cool postshock WR wind. The veloc- 
ity and density shear across the CD causes KH instabilities to 
develop and more- so in the trailing arm of the WCR due to the 
slightly higher shock obliquity close to the apex. As the stars 
move toward periastron the influence of orbital motion on the 
WCR becomes gradually more noticeable and the degree of cur- 
vature to the WCR arms increases. The curvature of the WCR 
reaches a maximum at periastron corresponding to the peak in 
the orbital velocities. Interestingly, as the stars approach perias- 
tron the postshock WR wind cools more rapidly, and the den- 
sity contrast across the CD increases. This, combined with the 
oblique angle at which the WR wind collides against the trail- 
ing arm of the WCR, results in perturbations to the CD growing 
to such an extent that filamentary structure is formed by com- 
pressions/rarefactions. We note that a similar feature wa s ob- 
served in simulations of rj Carinae by Parki n et all (1201 lb , but 
wa s not found in simulations of shorter period OB star systems 
bv lPittardl (120091). or in th e circular orbit simulations presented 
by Ivan Marie et al.l (1201 ll) . suggesting that it may only be ac- 
tive when the orbital eccentricity is relatively high, and hence 
the trailing arm of the WCR is considerably curved around pe- 
riastron. Although some of the postshock O star's wind remains 
quasi-adiabatic as it flows off the grid, the region close to the 
contact discontinuity cools more effectively. This is in part due to 
the mixing introduced by KH instabilities, but also results from 
some nume rical heat conduction b etween the hot and cold post- 
shock gas dParkin & Pittardll2Q10h . 

Following periastron, the WCR is wound around the stars by 
the influence of their orbital motion, and the cavity in the wake 
of the O star resembles an archimedian spiral (0 = 1.1 in Fi g.|4]) 
remin iscent of the so-called "pin wheel" nebula (e.g. lTuthilTe t al. 
2008). Examining the temperature snapshot at (p = 1.1 in Fig. [4] 
shows that postshock O star wind in the arms of the WCR cools 
to T - 10 4 K before exiting the grid. This gas was shocked at an 
earlier phase when xo was lower, however, xo was sufficiently 
large for the stars to have rotated in their orbits before cooling 
becomes apparent in the downstream gas. 

As the stars continue in their orbit the spiralling WCR grad- 
ually drifts outwards and exits the grid - a full rotation does not 
fit insi de the simulation domain. As shown by iParkin & Pit tard 
( 2008), the coils of the WCR will extend out to large distances 
(~ 1000 au). Heading towards apastron the orbital velocities de- 
crease and the impact of orbital motion on the WCR gradually 
fades, returning once again to an approximately axisymmetric 
shape. 

3.3. Model B Dynamics 

3.3.1. Large-scale phenomenon 

In model B the stellar winds are radiatively driven, and therefore 
consideration is given to their acceleration regions. Similar to 
model A, at apastron in model B the postshock winds are quasi- 



adiabatic and the effects of orbital motion are relatively minor. 
However, even at apastron the radiation fields of the stars intro- 
duce some noticeable differences. For instance, there is a small 
amount of radiative inhibition of the preshock WR wind by the 
O star which causes the distance between the WR and its respec- 
tive shock to be slightly smaller in model B compared to model 
A (Fig.0. Also, due to the lower off-axis velocities, the opening 
angle of the shocks are larger in model B compared to model A. 

Proceeding towards periastron, substantial differences be- 
tween models A and B become apparent. This is essentially 
tied to the increasing importance of cooling in the postshock 
gas as the WCR moves deeper into the wind acceleration re- 
gions (§ I3.ll) . For instance, between ^ 0.8 - 0.9 there is a 
sharp increase in the importance of cooling in the postshock 
O star's wind and consequently thin shell instabilities become 
more apparent at the apex of the WCR. The increasing level of 
radiative inhibition, which occurs as the stellar separation con- 
tracts, assists the disruption of the WCR by further reducing the 
O star's preshock wind velocity. Interestingly, the cooling pa- 
rameters calculated from model B agree very well with the es- 
timates in § I3.ll based on approximating the wind acceleration 
using /3-velocity laws (Fig. [2]). There is, however, a spike at 
cp ^ 0.93 corresponding to instabilities perturbing the position 
of the WCR and thus the acquired preshock velocities. 

As the stellar separation contracts, the WCR moves closer 
to the O star and its postshock wind begins to cool effectively, 
forming a thin, dense layer. When both win ds cool effecti vely 
postshock, thin-shell instabilities at the CD ( Vishniacl ll983l) are 
permitted to grow in amplitude. Essentially, this occurs because 
the restoring force provided by thermal pressure in the postshock 
gas becomes ins ufficient to preve nt bending modes from evolv- 
ing non-linearly (IVishniaclll994l) . This contrasts with model A, 
where the postshock O star's wind is quasi-adiabatic throughout 
the orbit, and the WCR is supported against the growth of NTSIs 
because the high thermal pressure of the postshock O star wind 
acts like a cushion and suppresses the growth of ripples in the 
CD. 

The onset of effective radiative cooling for both winds in 
model B also considerably changes the downstream flow in the 
WCR arms, making it far more structured/clumpy than in model 
A (see cp = 1.1 snapshots in Figs. |4] and O. During the tran- 
sition from quasi-adiabatic to radiative postshock gas the com- 
bination of cooling and compression by KH instabilities forms 
clumps in the downstream flow (see cp = 0.90 snapshot in Fig. [6]). 
Moreover, the vigorous action of the NTSI breaks-up the thin- 
shell of postshock gas as it flows away from the apex of the 
WCR. Interestingly, although the NTS I was found to develo p 
in the simulations of iPittardl (l2009b and Ivan Marie et al.l (120111) . 
such excessive fragmentation was not observed, perhaps due to a 
lower resolution of the postshock flov\0. Before (after) periastron 
the trailing (leading) arm of the WCR is more structured than the 
leading (trailing) arm due to the larger shock obliquity and thus 
more effective growth of KH instabilities (see the cp = 0.97 and 
1.05 snapshots in Figs. [6] and [7]). 

As the stars proceed towards apastron the continuing injec- 
tion of thermal pressure at the apex gradually restores the sta- 
bility of the downstream shocks and smoothes out the clumpy 
structure in WCR arms. As the stars proceed towards apastron 
the postshock winds become quasi-adiabatic again. The rela- 
tively rapid destruction of clumps in the WCR following peri- 



5 As demonstrated by IParkin & Pittardl (|2010), the degree of frag- 
mentation of the thin shell increases, and the size of the resulting clumps 
decreases, with the increasing resolution of the postshock flow. 
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Fig. 4. Snapshots of the gas density (left column) and temperature (right column) in the orbital (x - y) plane from model A (instan- 
taneously accelerated winds) at cp = 0.5 (upper panels), 1.0 (middle panels), and 1.1 (lower panels). The orbital motion of the stars 
is calculated in the centre of mass frame. At apastron (0 = 0.5) the WR star is to the right, and the O star is to the left, of the image 
centre. The motion of the stars proceeds in an anti-clockwise direction. All plots show a region of ±1.2 x 10 14 cm - large axis tick 
marks correspond to a distance of 5 x 10 13 cm. 



astron contrasts with the simulation s of shorter period OB star 
systems presented by iPittardl (|2009), in which clumps formed 
close to periastron survived until apastron. 



3.3.2. The collapse/recovery of the wind-wind collision region 

The growth of NTSIs corrugates the WCR apex, increasing the 
shock obliquity, thereby reducing the efficiency with which wind 



kinetic energy is thermalized. Consequently, postshock gas tem- 
peratures are further reduced and cooling becomes even more 
effective. The runaway disruption continues as periastron is ap- 
proached and the WCR apex oscillates wildly in the vicinity of 
the O star until a collision occurs at ^ 0.95. The preshock 
O star's wind, which suffers substantial radiative inhibition at 
phases close to periastron, has little room to accelerate and its 
ram pressure is unable to hold back the incoming WR wind, 



6 



E. R. Parkin and E. Gosset: WR 22: Hydrodynamical modelling and X-ray emission 




Fig. 5. Same as Fig.|4]except model B (radiatively driven winds) is shown. 



leading to a collaps^l of the WCR onto the O star at = 0.96. 
We note that contrary to the predictions of Fig. |3]radiative brak- 
ing does not prevent a collapse of the WCR because the region 
over which braking is predicted to occur is occupied by post- 
shock gas at T > 10 6 K which we assume to be mostly ion- 
ized and consequently to be subject to a negligible driving force. 



6 During a collapse of the WCR onto the O star, the WR's wind will 
enter deep into the O star's wind acceleration region, and may poten- 
tially collide against the O star's photosphere. We do not expect that 
this will lead to significant accretion and mass transfer as the WR's 
wind has too great a kinetic energy to become bound to the O star (c.f. 

§ED. 



How ever, the strength of radiative braking is dependent on th e 
ICAKl parameters adopted, i.e. the k and a (iGayley et al.lll997l) . 
When calculating the line force appl ied to t he WR's wind by the 
O star's radiation we use the O star's ICAKl parameters (ko = 0.3 
and ao = 0.52). Hence, if the WR's radiation-wind coupling 
was adopted (&wr = 0.42 and #wr = 0.47) the decelerative line 
force would be greater and radiative braking may be more effec- 
tive than in our calculations. 

Interestingly, as the stars reach cp = 1.0 the apex of the WCR 
is essentially stabilised against NTSIs by the continuing collapse 
(i.e. the ram pressure of the WR's wind pins the WCR apex to the 
O star preventing oscillations) and only relatively small ampli- 
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tude KH instabilities can be seen in the trailing arm of the WCR 
(Fig. [(J). Despite the O star's wind being overwhelmed between 
the stars it continues to drive a wind on the side facing away from 
the WR and, as in model A, carves a tenuous cavity through the 
dense WR wind (Fig. 0). As the wind driven from the far side 
of the O star shocks against the leading arm of the WCR (which 
becomes wrapped around the stars) it reheats the postshock gas, 
and consequently the leading arm of the WCR is hotter than the 
trailing arm (see the temperature plots at = 1.0 and 1.1 in 
Figs.@]and[5]). 

Following periastron, as the separation of the stars expands, 
the preshock ram pressure of the WR's wind - which pins the 
shocks onto the O star - decreases. Therefore, when it weak- 
ens sufficiently, the WCR is permitted to move away from the 
O star. Simultaneously, as the stars recede, radiative inhibition 
decreases and the acceleration of the O star's wind increases. 
The O star shock is buried so deeply in its wind acceleration 
region that a small difference in its position considerably alters 
the preshock velocity. Hence, as soon as the O star's wind is 
given some room to accelerate the acquired preshock ram pres- 
sure greatly increases, helping the WCR to separate from the O 
star. Subsequently, the higher preshock velocities lead to higher 
postshock gas temperatures/pressures, which increase the stabil- 
ity of the WCR apex against NTSIs. The snapshots in Fig. [7] 
show the gradual transition from collapse to recovery. Notice 
that the recovery is not instantaneous and, as with the collapse, 
NTSIs vigorously disrupt the WCR and sporadic collisions of 
the WCR against the O star occur (e.g. = 1.08 in Fig. [7]). 
However, as the amplitude of the oscillations increases, the O 
star's wind is allowed to accelerate to higher velocities before it 
shocks and, as such, the postshock gas pressure sees a commen- 
surate increase. Following (p = 1.10 the WCR maintains a clear 
separation from the O star. 

On the basis of the wind-wind momentum balance alone, 
the separation of the WCR from the O star is inevitable (Figj2]). 
However, by initiating the winds in model B wi thin a radial dis- 
tance of ~ 1.15 R* (where 0.15 R* > 3 cells - lPittardll2Q09b we 
are essentially enforcing a lower limit to the wind ram pressure 
(i.e. the O star's wind velocity never reaches zero) and influenc- 
ing the time at which the shocks d etach from t he O star. For 
instance, in models of Iota Orionis, Pittard (1998) found that the 
orbital phase at which the WCR recovered was tied to the size of 
the annulus used to initiate the stellar wind with a smaller annu- 
lus leading to a more prolonged collapse. In reality, the incoming 
wind may interfere with the initiation of the wind and prolong a 
collapse/disruption of the WCR. Hence, although the exact time 
at which the WCR separates from the O star may be resolution 
dependent, its occurrence in model B is qualitatively correct. 

We note that numerical heat conduction between hot ten- 
uous gas and cool de nse gas will occur in the simulations 
dParkin & Pittardl 12010b . In model B this will affect the tem- 
perature of postshock gas, and hence the rate at which it cools. 
For instance, during the collapse of the WCR the postshock WR 
wind will be caused to cool more rapidly and the O star wind 
will be artificially heated, potentially to X-ray emitting tempera- 
tures. This effect is also present in model A, but to a lesser extent 
due to the lower contrasts between adjacent postshock gas (i.e. 
there are cool clumps in the postshock gas in model B but not 
in model A). However, the conclusions drawn from this work 
are not affected by the undesirable influence of numerical heat 
conduction. 
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Fig. 6. Sequence of snapshots of gas density (left column) and 
temperature (right column) in the orbital (x-y) plane illustrating 
the collapse of the WCR onto the O star. From top to bottom: 
= 0.90, 0.92, 0.95, 0.97, and 1.00. All plots show a region of 
±3 x 10 13 cm. 
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Phase = 1 .08 




Fig. 7. Same as Fig |6] except showing the recovery of the WCR 
from collapse. From top to bottom: = 1.05, 1.06, 1.08, and 
1.10. All plots show a region of +3 x 10 13 cm. 



3.4. X-ray emission 

The collision of the winds in WR22 generates hot (T > 10 7 K) 
plasma which emits radiation at X-ray wavelengths, providing 
a direct probe of the postshock winds, and an indirect probe of 
the preshock winds. Howe ver, a recent analysis of XMM-Newton 
observations of WR22 by lG09l revealed a number of discrepan- 
cies between predictions from a wind- wind collision model and 
observations. To determine whether the contrasting flow dynam- 
ics from models A and B can shed any light on these issues we 
have performed X-ray radiative transfer calculations on the sim- 
ulation output, details of which can be found in § 12.31 We define 
our line of sight geometry as follows: the inclination angle, i, is 
measured against the z-axis (z = 0° would view the system from 
directly above the orbital plane), and the angle 6 is measured 
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Fig. 8. X-ray lightcurves calculated from models A and B in the 
0.5-1 (upper panel), 1-2 (middle panel), and 2-10 keV (lower 
panel) energy bands. The values derived by lG09l from spectral 
fits to XMM-Newton spectra are plotted for comparison - note 
that they have been rescaled (upwards) by a factor of 200 for 
illustrative purposes. Note the difference in scale between the 
plots. 



against the negative x-axis (O star in front at apastron) such that 
6 increases in the prograde direction (0 = 90° would align the 
line of sight with the negative y-axis). We adopt viewing angles 
of i = 85° and = 0°, in agreement wit h the system orienta- 
tions d etermined by lRauw et al.1 (Il996l) and lSchweickhardt et al.l 
(119991) . 
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3.4.1 . Lightcurves 

The 0.5-1, 1-2, and 2-10 keV X-ray lightcurves calculated from 
models A and B are shown in Fig. [8] Although the separate con- 
tributions are not shown, the X-ray emission is dominated by the 
postshock WR wind due to its higher density and mean atomic 
weight. Examining the lightcurves for model A one sees a rea- 
sonably constant flux between 0.35 - 0.65 corresponding to 
the relatively slow rate of change of the binary separation whilst 
the stars are close to apastron. As the stellar separation contracts, 
the preshock, and thus postshock, gas density increases resulting 
in a rise in the observed flux. However, approaching periastron, 
lines of sight to the X-ray emitting plasma begin to pass through 
the WR's wind and the abrupt increase in absorption causes the 
observed flux to decline. The turn-over occurs at an earlier phase 
for X-rays with E < 2 keV due to the higher susceptibility to 
absorption compared to those at E > 2 keV. The decline in ob- 
served flux continues until a minimum is reached at periastron 
passage thence attenuation reaches a maximum. The spatial ex- 
tent of the X-ray emission region means that an eclipse by the 
WR star will have only a small impact, hence the dip in ob- 
served flux around periastro n mainly occurs due to absorption 
by the unshocked WR wind dParkin & Pittardll2008h . 

The model B lightcurves differ from those of model A in 
a number of ways. Firstly, the observed flux at apastron in 
model B is higher than in model A due to the distance be- 
tween the WR and its respective shock being slightly smaller 
in the former (see § 13.31) . Secondly, the decline in flux begins 
at the same orbital phase (0 ^ 0.8) in the 0.5-1, 1-2, and 2- 
10 keV lightcurves. Interestingly, in model B the start of the de- 
cline in the 2-10 keV flux corresponds to the initiation of the 
WCR disruption, unlike in model A where no disruption occurs. 
Similarly, the egress in the observed flux following periastron 
is tied to the recovery of the WCR. The highly unstable WCR 
introduces flare-like features into the lightcurves. For example, 
during the rapid decline in the observed flux (0 ^ 0.8 - 1.0) 
there is a sudden spike at = 0.94 caused by the WCR os- 
cillating away from the O star (Fig. [6]). It is interesting to note 
that flare-l ike features are observed in the X-ray lightcurve of 
rj Car inae (iCorcor an et al. 2001; Corcoran 20051 ICorcoran et aD 
12010b . Although [Moffat & Corcoran! (120091) dismissed instabili- 
ties in the WCR as the mechanism responsible for the flare-like 
features, the model B lightcurve provides some evidence to the 
contrary. 

The most apparent difference between models A and B is the 
drastic reduction in the 2-10 keV flux caused by the disruption, 
and subsequent collapse onto the O star, of the WCR in the lat- 
ter. Consequently, at periastron the difference between the two 
models reaches a maximum with the 2-10 keV flux from model 
B being a factor of ^ 1 8 lower than from model A. This occurs 
because during the collapse the fraction of the WR's wind being 
shocked is approximately the fractional solid angle subtended by 
the O star, which is much lower than occurs when a stable WCR 
exists. 

For comparison, we also plot the X- ray fl uxes derived from 
XMM-Newton observations of WR22 by G09 in Fig. [8] Clearly, 
the models o ver-p redict the G09 values in all three energy bands 
(note that the lG09ll values have been rescaled in Fig. [8]). At apas- 
tron, when a stable WCR exists in both models, the 0.5-1, 1-2, 
and 2-10 keV fluxes are over-estimated by factors of 220, 230, 
and 170, respectively. The major concern is the over-estimate in 
the 2-10 keV flux, as X-rays in this energy band are less suscep- 
tible to absorption and, therefore, the discordance between the 
models and observations cannot be attributed to insufficient at- 
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Fig. 9. Emission weighted column density ( EWC) calculated 
from models A and B. The values derived by G09 from spec- 
tral fits to XMM-Newton spectra are also plotted (see § 13.4.2b . 
with the formal error bars attained from the spectral fitting . The 
plotted values do not include the column density due to the in- 



terstellar medium (2.5 x 10 21 cm z ). 




Fig. 10. Emission weighted column density (EWC) as a function 
of energy calculated from models A and B at = 0.994. The 
column density due to the interstellar medium (2.5 x 10 21 cm -2 ) 
is additional to the plotted values. 



tenuation. Hence, to reduce the observed 2-10 keV flux requires 
a reduction in the intrinsic flux, i.e. an alteration to the WCR. 
Encouragingly, the disruption/collapse of the WCR closes the 
gap considerably; at = 0.994 model B only over-estimates the 
G09 values by factors of ^ 7.8, 4.7, and 6.8 in the 0.5-1, 1-2, 
and 2-10 keV energy bands, respectively. The comparison of the 
model lightcurves and the G09 values would therefore suggest 
that the WCR must be significantly disrupted/collapsed onto the 
O star throughout the entire orbit. We consider this point further, 
along with revisions to parameters to facilitate it, in § [ 
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3.4.2. Column density 

The emission-weighted column densities (EWCsfl calculated 
from models A and B are plotted in Fig. [9] A minimum in the 
column density occurs for both models at = 0.5 corresponding 
to the stars at apastron and the lines of sight to the X-ray emit- 
ting plasma passing through the O star wind. The column density 
then rises gradually between = 0.5 - 0.8 as the WCR, and thus 
lines of sight to the X-ray emitting plasma, move deeper into the 
O star wind. At = 0.8 the column density begins to rise rapidly 
due to the WCR rotating and some lines of sight now viewing 
the X-rays through the WR wind. The column density peaks at 
periastron when the WCR is at its deepest point of immersion 
in the dense WR wind, following which there is a decrease as 
the stars recede and lines of sight to the X-ray emitting plasma 
traverse progressively lower density wind material. Comparing 
the present model column densities to those computed with the 
model of lG09l a good agreement is observed. 

The general morphology of the column densities calculated 
from models A and B are largely similar for most of the orbit, 
however, there are some noticeable differences. The column den- 
sity is higher for model B compared to model A at = 0.5 due 
to the wind density being slightly higher for the former in the 
unshocked O star wind immediately adjacent to the WCR apex 
(where the bulk of the X-ray emission originates from) - this is 
a numerical artifact of the different ways by which the winds 
are driven in the two simulations. The step in column density in 
model B which extends to ^ 1.08 corresponds to lines of sight 
to the X-ray emitting plasma passing through cool, dense post- 
shock gas (see Fig. [7]). Interestin gly, a similar feature is observed 
in si mulations of rj Carinae (Park in et~aD l2009al iParkin et aD 
[2oTTh . As the O star wind recovers, the postshock gas temper- 
ature suddenly increases and cooling becomes less important, 
thereby reducing the postshock gas density. Simultaneously, the 
rotation of the arms of the WCR results in fewer lines of sight 
to the X-ray emitting plasma being obscured by dense postshock 
gas. 

For comparison, two sets of the column densities derived by 
G09 are also plotted in Fig.[9j the lower values correspond to fits 
performed with a single column density to account for all ab- 
sorption (G09 - soft Nu - see table 7 of lG09l) . whereas the higher 
values correspond to the harder component in a two-temperature 
fit w ith se parate column densities (G09 - hard Nu - see table 
6 of lG09h . We note that the "s oft" c olumn densities are rep- 
resentative of values derived by lG09l to the lower temperature 
plasma component in two-temperature, two-column density fits. 
The observed column density to the hotter emission component 
is higher than to the softer plasma component (although these 
values could be lower - see lG09l) . the physical interpretation of 
which would be that the hotter plasma resides closer to the WCR 
apex and hence lines of sight pass through denser wind mate- 
rial. In contrast, the lower column density accrued along lines of 
sight to the cooler temperature plasma suggests that this emis- 
sion comes from a more extended region which is largely viewed 
through less dense wind material. Unsurprisingly, the column 
densities calculated from models A and B lie roughly between 
the two sets of observed column densities as they sample a range 
of values encountered when viewing a spatially extended, energy 
dependent emission region. Fig.[l0]shows the column density as 
a function of energy at = 0.994 - the column density at 10 keV 

7 The EWC is calculated as N E wc = ^Mi^intx/^intx, where Nu and 
L int x are the column density and 0.5-10 keV intrinsic luminosity from a 
given line of sight, a nd the summation is over all sight lin es (pixels) in 
the X-ray image fe.g. lParkin & Pittard 2008; Parkin et al. 200§B). 



is an order of magnitude higher than at 0.5 keV. This highlights a 
difficulty when applying simple spectral fitting models to com- 
plicated data, namely th at observationally derived values may 
be so mewhat misleading (lAntokhin e t al. 2004t lPittard & Parkinl 
l2010h . For exam ple, b ased on an over-prediction of the observed 
column density, lG09l stated that a more extended emission re- 
gion was required to reduce the model column density and attain 
a better agreement. Yet, as demonstrated by Fig. [9] either an ex- 
tended (model A) or relatively small (model B) spatial extent to 
the emission region at = 0.994 both produce column densities 
which lie between the observed values derived by lG09l for the 
"soft" and "hard" hot plasma components. 

Comparing the column densities from models A and B one 
can see that the E > 2 keV values are higher in the latter, which 
results from the smaller spatial extent of the emission region at 
this orbital phase in model B compared to model A (see Figs. [4] 
and0. 

3.4.3. X-ray spectra 

The observed X-ray spectrum from the WCR is the product of 
viewing a spatially extended region which emits at a range of 
energies through an intervening medium (i.e. the line of sight 
absorption), which may itself have a complicated spatial depen- 
dence. Therefore, despite the models over-predicting the inte- 
grated 0.5-10 keV flux, a comparison of the shape of the model 
spectra against observations can provide useful insight. The X- 
ray spectra calculated from models A and B at = 0.63, 0.903, 
and 0.994 are plotted in Fig. [TT1 At 6 = 0.63, a stable WCR is 
established in models A and B, and consequently the X-ray spec- 
tra calculated at this phase appear largely similar (Fig. [TT|). Both 
mod els provide a good fit to the = 0.582/0.680 spectrum from 
G09, albeit with a significant overestimate in the flux level (note 
that the model spectra have been rescaled to have an equivalent 
10 keV flux). As one progresses closer to periastron, the agree- 
ment between the shape of the model spectra and observations 
worsens; at = 0.903 the models underpredict the E ^ 1 keV 
flux, whereas at = 0.994 the observed spectrum at E < 8 keV 
is poorly matched by either mode@. In particular, there appears 
to be too much absorption at E < 4 keV close to periastron, 
consistent with the findings of lG09l 

The lack of agreement between the shape of the model A and 
B spectra and the lG09l spectrum at = 0.994 may be indicating 
that our adopted viewing angles are incorrect and that some ad- 
justments are warranted. To explore this possibility tests have 
been performed with viewing angles in the range i = 70° to 90° 
and = -10° to 40°. In essence, the degree of absorption 
can be reduced by adopting viewing angles which peer through 
less dense WR wind material. However, although the amount of 
E < 4 keV emission is slightly increased, a significant discrep- 
ancy still remains. To illustrate this point the model A and B 
spectra calculated using i = 70° and = 20° at = 0.994 are 
also plotted in Fig. [TU the improvement in the match is small 
and the observed E = 2 keV flux (in the normalized spectrum) 
is underestimated by over two orders of magnitude. 

The deficit in E < 4 keV flux (in the normalized spectrum) 
at = 0.994 could be partially remedied by considering the ion- 



8 The temperatu re of the hotter plasma component is undefined in 
the spectral fits of G09 at = 0.994. Therefore, to facilitate a compar- 
ison between the models and observations, we take kT 2 = 4 keV and 
Norm<i = 1.1 x 10~ 4 cm -5 which leads to an observed 0.5-10 keV flux of 
^ 1.2 X 10~ 13 erg s _1 cm -2 , in agreement with, but slightly lower than, 
table 12 oflG09l. 
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The X-ray spectra also provide a direct probe of postshock 
gas temperature, and therefore the preshock wind velocity. In 
particular, the slope of the spectrum at E > 5 keV is a useful 
indicator of the validity of an adopted wind velocity as X-rays 
at this energy are less susceptible to absorption. Examining the 
spectra at = 0.582/0680, 0.903, and 0.994 one sees that model 
A provides a slightly better match to the slope of the spectrum 
at E > 5 keV than model B. Considering that at = 0.994 in 
model A the WR wind has a preshock velocity of 1785 km s" 1 , 
whereas in model B radiative inhibition causes a slight reduction 
down to 1660 km s" 1 , this could be indicating one of two things. 
Firstly, the terminal velocity of the WR wind is in fact higher, 
and that at = 0.994 radiative inhibition reduces this velocity 
to a preshock value of 1785 km s" 1 . However, this seems un- 
likely as the termina l wind speed of the WR wind is known wit h 
some confidence (Crowthe r et al.l [l9 95a; Hamann et al. 2006). 
Alternatively, radiative inhibition for WR win ds is not as ef- 
fective as predicted by f ormulations based on ICAKl theory (i.e. 
IStevens & Pollocklll994l) . Detailed non-LTE models of WR1 1 1 
support this idea as they find the lin e force to be insensitive t o 
the gradient in the wind velocity (Grafener & Ha mannll2005l) . 
Furthermore, the ionizing radiation from the O star may dramat- 
ically increase the line force in the incoming WR wind, hence 
considerably enhancing the strength of radiative braking and af- 
fecting th e opening angle of the WC R and the resulting X-ray 
emission (Grafener & Hamann 20071). 
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Fig. 11. 0.5-10 keV X-ray spectra calculated from models A 
and B at = 0.630 (upper panel), 0.903 (middle panel) and 
0.994 (lower panel). Model spectra calculated at cp = 0.994 us- 
ing viewing angles of i = 70° and 6 = 20° are also shown in 
the lower panel. The best-fit model spectra from lG09l are plot- 
ted for comparison (see table 7 of G09). Note that the average 
cp = 0.582/0.680 spectrum of G09 is plotted in the upper panel. 
For illustrative purposes, the model A and B spectra have been 
rescaled by the factors noted in the plots to have an equivalent 
10 keV flux. Note the difference in scale between the plots. 



ization state of the inner WR wind. For instance, if the gas is 
highly ionized it will be less opaque. However, opposite con- 
straints could come from the analysis of the eclipses at optical 
wavelengths. In any case, an in-depth examination of this effect 
is beyond the scope of the current work. 



3.5. A revised model for WR 22 

The comparison between the model results and observations in 
the previous sections highlights a number of discrepancies, the 
most severe being a substantial overestimate of the X-ray flux 
across the majority of the orbit. The major concern is the over- 
estimate of the 2-10 keV X-ray flux as, compared to X-rays at 
E < 2 keV, this emission is less affected by absorption, and 
therefore a large reduction in the intrinsic emission is required. 
In this section we consider revisions to model parameters which 
may improve the agreement with observations. 

The collapse of the WCR in model B considerably improves 
the agreement with the observed X-ray flux. Essentially, this oc- 
curs because the fraction of the WR's wind being shocked is 
vastly reducecQ. As a result, a sustained collapse/disruption of 
the WCR could provide a better agreement between the model 
and observations. Such a situation would occur if the ram pres- 
sure of the WR's wind massively overwhelms that of the O star's, 
which could be engineered in a number of ways. Firstly, one 
could increase the mass-loss rate of the WR's wind. Noting that 
the WR mass-loss rate adopted in this work is relatively low in 
comparison to the ^ 4x 10~ 5 M n yr" 1 and 6.3x!0~ 5 M vr" 1 de - 
termined by ICrowther et al.l (Il995al) and lHamann et al.l (120061) . 
respectively, this would be a reasonable alteration, but at the bor- 
der of the possibilities authorized by the eclipses in the visible 
domain. Alternatively, or simultaneously, one could reduce the 
mass-loss rate of the O star. Considering the uncertainty in the 
parameters of the O star, a factor of ~ 2 - 3 lower mass-loss rate 
is within reason. Repeating the calculations in § 13.11 in which 
the winds are approximated by ft- velocity laws, we find that if 
the ratio of the WR and O star mass-loss rates is increased from 
^ 57 (TableO to ^ 500, there would be no stable wind- wind mo- 



9 By initiating the winds within a radius of ~ 1.15 R* we are artifi- 
cially increasing the fraction of the WR wind being shocked by essen- 
tially augmenting the size of the O star. However, the additional X-ray 
flux introduced has a negligible effect on the results. 
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mentum balance at any orbital phase. Considering that a larger 
WR mass-loss rate may increase the amount of absorption at pe- 
riastron, and therefore potentially worsen the fit of the model 
spectrum to observation, the most suitable option would appear 
to be a reduction of both the WR and O star mass-loss rates, with 
the ratio of the two being ^ 500. 

If the WR's wind ram pressure significantly overwhelms 
the O star's then r adiative braking could play a significant role 
(iGayley et al.ll 19971) . Fig. [3] shows that, using the stellar param- 
eters in Table |2 the O star should fully brake the incoming WR 
wind at apastron, but may only partially brake it at periastron. 
This would remain the case if the WR's mass-loss rate is in- 
creased as proposed above. Consequently, the relatively small 
increase in the observed 2-10 keV flux from apastron to perias- 
tron may be due to the WCR transitioning from a ram balance 
(supported by radiative braking) to a collision against the O star. 
Curiously, the good agreement between the model spectra and 
observations at = 0.582/0.680 could be indicating that the 
WCR apex may not be collapsed onto the O star at this phase, 
perhaps supported by effective radiative braking at phases close 
to apastron. However, if a WCR were established there would 
either need to be a significantly lower fraction of the WR's wind 
being shocked, or much lower density postshock gas, to avoid 
the resulting X-ray flux exceeding observations. 

There is also the question of whether the agreement be- 
tween the models and observations could be improved if wind- 
clumping is considered. Examining Fig. [8] at cp ^ 0.58 the 
models overestimate the observed 2-10 keV flux by a fac- 
tor of ^ 170 which, as Lx oc M 2 , implies a reduction in 
mass-loss rates by a factor of ~ 13. This greatly exceeds 
the current estimates of reduction factors of 2- 5 for global 
wind mass-loss rates to account for cl umping (iBouret et al.l 
2003; Repolu st et al.l|2004l: iMarkoya et al.ll2004l: iFullerton et al.1 
2006; Crowth edl2007l;lMoffatll2008l:lPuls et al.ll2008l) . Therefore, 
clumping alone cannot account for the discrepancy between the 
models and observations. 

In summary, the model X-ray flux will likely agree better 
with observations if the WR's wind ram pressure overwhelms 
the O star's. This could lead to a collapse at all phases which 
would drastically reduce the observed X-ray flux (e.g. model B 
at phases close to periastron in Fig. [8}. There is, however, the 
intriguing possibility of a WCR established at apastron by radia- 
tive braking transitioning to an instability-driven disruption, or 
collapse of the WCR onto the O star, at periastron. For instance, 
the goodness of fit of the spectral shape at apastron argues for 
a stable WCR which, in the absence of a normal ram pressure 
balance, could be established by radiative braking. Yet, although 
a collapse provides better agreement in terms of the integrated 
flux, the shape of the model B spectrum at cp = 0.994 does not 
agree very well with observations as there is too much absorp- 
tion at E < 2 keV. As already mentioned, this may be indicat- 
ing that the WR's m ass-lo ss rate requires a downward revision 
as also suggested by lG09l . However, a minor improvement may 
also be attained by considering the high ionization state, and thus 
lower opacity, of the obscuring inner WR wind. Furthermore, 
if radiative braking is effective at periastron then the increased 
opening angle of the WCR may result in a sufficiently extended 
emission region such that some E < 2 keV X-rays pass through 
less dense WR wind material. Minor reductions in the model flux 
could also be achieved by adopting a smaller O star radius and/or 
reducing the global mass-loss rates in-line with the range of un- 
certainty permitted by current estimates of wind-clumping. The 
former could reduce the model flux during a WCR collapse by 
reducing the fraction of the WR's wind being shocked, whereas 



the latter would lower the emission measure of the postshock gas 
at all orbital phases. 

4. Discussion 

The simulations presented in this work, although detailed and 
insightful, represent an approximation to reality. In this section 
we consider the potential importance of physics and/or processes 
which may not be captured by our simulations, and an inherent, 
although minor in this case, drawback of a numerical approach. 
These estimates are by no means a substitute for further detailed 
investigations, and are merely meant to serve as an indicator that 
there is no obvious candidate responsible for the discrepancy be- 
tween the observed and model spectra at cp = 0.994 (Tig.fTTl). 

4.1. Accretion 

During periastron passage in model B, a stable WCR is not es- 
tablished due to the overwhelming dominance of the WR's wind 
over the O star's wind. In § !3.3.2l we describe this as a collapse 
of the WCR onto the O star, and it may therefore seem reason- 
able to consider the possibility that the O star may accrete some 
of the WR's wind. 

Capture of accreted mass will occur within the Hoyle- 
Lyttleton radius dFrank. King & Raindl2002h . 

2GM 



? ? ' 

V M + V™ 

*0 WR 



(5) 



where v*o and vwr are the orbital velocity of the O star and wind 
velocity of the WR wind, respectively. Inserting values appropri- 
ate for WR22 at ^ 1.0, and assuming the WR wind to be at 
terminal velocity, we find r acc 3 R , noticeably smaller than 
the adopted O star radius (Table 13. As r acc gives an approx- 
imate value for the radius at which an incident flow becomes 
gravitationally bound, this result suggests that, despite the close 
proximity of the WR's wind to the O star, significant accretion is 
unlikely. One may argue that if the incoming WR wind is suffi- 
ciently radiatively braked it may lead to some capture. However, 
setting r acc = Ro and rearranging Eq. |5]to determine the required 
O star velocity for accretion to occur it becomes apparent that 
vwr must vanish. As Fig. [3] shows, for our adopted radiation- 
wind coupling this does not appear to be the cas<o» With these 
points considered, we do not anticipate that accretion will oc- 
cur in WR22 during a collapse of the WCR onto the O star at 
periastron passage. 

4.2. Thermal electron heat conduction 

In reality, some thermal conduction of heat across the con- 
tact discontinuity may occur, yet its magnitude is depen- 
dent on the magnetic field s trength and orienta tion (see e.g. 
Mvasnikov & Zhekov] 119981: [Orlando et al. 2008). If thermal 
electron heat conduction were to be effective, the postshock gas 
temperature c ould be reduced, softening th e intrinsic X-ray spec- 
trum (see e.g. lZhekov & Mvas nikov 2003). An indication of the 
importance of thermal conduction in the WCR shocks can be 
gained from an exa mination of the thermal conduction timescale 
( Orl ando et aDl20Q5h . 



^cond 



ps 



2( r -i)^(r ps )(r ps // 2 ) 



s, 



(6) 



10 However, if the strength of the decelerative force of the O star's ra- 
diation is far greater (c.f. iGrafener & H amann 2007) then the incoming 
WR wind may be braked more effectively. 
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where P ps and T ps are the postshock gas pressure and temper- 
ature, respectively, and / is a characteristic length scale (taken 
here to be the width of the shock at the apex of the WCR in the 
simulations). Note that we assume a single temperature for the 
postshock electrons a nd io ns which should pr ovide an accurate 
approximation (c.f. § 14.41) . The conductivity (iBorkowski et al.l 
11990b . 

T 5/2 

k(T vs ) = 5.6 x 10" 7 erg cm -1 s" 1 K" 1 , (7) 

v p J 1 + (£tan<9) 2 

where f is the ratio of the number density between the cool and 
hot gas and 6 is the angle between the magnetic field direction 
and the shock normal. Evaluating Eq. [6] for the case of 6 = 0° 
(i.e. magnetic field lines aligned parallel to the shock normal), 
and taking the value of / from model B, we find £ con d - 4.7 days, 
which is notably longer than the timescale for postshock gas to 
advect away from the apex of the WCR at phases close perias- 
tron, r gas = J/vwr - 0.7 days, where d is the binary separation. 
Examining the alternate extreme, where the magnetic field lines 
are approximately perpendicular to the shock normal, we find 
^cond - 40 and 9900 days for thermal conduction between the 
postshock WR wind and the preshock WR wind and the post- 
shock O star's wind, respectively. As a result, based on the above 
estimates, we do not anticipate that thermal conduction will be 
largely important for explaining the discrepancy between the ob- 
served and model spectra at = 0.994. 

4.3. Particle acceleration 

A fraction of the available particles could be accelerated at 
the shocks to re lativistic energies (for a recent review see 
iDe Beckerl 120071) . and in so doing the fraction of thermally 
emitting particles would be reduced. The efficiency of par- 
ticle acceleration in the wind-wind collision shocks of mas- 
sive star binary systems is unclear, with estimates ranging 
from < a percent (Eichler & UsovlH993t Dougherty et al. 2003; 
Pitta rd & Doughertvl l2006h up to ~ 50 percent dPittard et al.l 
2006^ If particle acceleration were to be efficient, shock mod- 
ification may occur whereby the diffusion of non-thermal ions 
upstream of the subshock exert a back-pressure on the preshock 
flow, causi ng the gas velocity to de crease prior to the sub- 
shock (see iPittard & Doughertyll2006l and references there-in). 
Consequently, the postshock temperature would be reduced and 
the observed spectrum softened. Such shock modification is, 
however, unlikely in WR22 at periastron due to the close prox- 
imity of the shocks to the O star, whereby inverse Compton 
losses will limit the growth of Lorentz factors for postshock 
ions. Furth ermore, WR22 has been cla ssified as a thermal ra- 
dio source (Doughe rty & Williamsll200(j) . Therefore, attributing 
the difference between model and observation at (p = 0.994 to 
non-thermal effects lacks any clear observational driver. 

4.4. Non-equilibrium ionization 

In the models presented in this work we assume that the 
postshock ions and electrons are in temperature equilibrium. 
Interestingly, studies of the massive star binary systems WR140 
and WR147 find that the observed soft X-ray continuum emis- 
sion can naturally be exp l ained by non-equili b rium ionizatio n 
(IZhekov & Skinned I200QI: iPollock .et al.l 120051: IZhekovl 120071) . 
namely the fact that temperature equilibration between post- 
shock ions and electrons via Coulomb collisions may occur at 
some distance downstream of the WCR apex. There are, how- 



ever, two important points which suggest that non-equilibrium 
ionization is not the missing ingredient in our models of WR22. 
The first is that in WR140 and WR147 the shocks are thought to 
be collisionless. In contrast, at periastron in WR22, the length- 
scale for Coulomb c ollisional dissipation, k-{ ^ 7 x 10 18 Vg/ft; 
(Pollock et al . 2005), where n x is the ion number density, is 
^ 4 x 10 9 cm, i.e. smaller than the finest simulation grid cell. 
The second point relates to the time scale for tem perature equili- 
bration between electrons and ions (Spitzer 1962), 

^ eq ~ 252^(lnA)- 1 -^s, (8) 

where \i is the mean molecular weight, Z is the metallicity, and 
In A is the Coulomb logarithm. Evaluating Eqj8] for postshock 
gas conditions in WR22 at periastron we find t Qq 0.2 days, 
which implies that the electrons and ions will rapidly equilibrate 
their temperatures close to the WCR apex. 

4.5. Numerical heat conduction 

Highly radiative shocks in colliding winds binary systems are 
difficult to model dMvasnikov et al. Il 19981: lAntokhin etaDl2004l: 
iParkin & Pittardl l2010h . This is mainly because of a significant 
difference in length scales; the cooling length for radiative post- 
shock gas may become many orders of magnitude smaller than 
the binary separation. Hence, sampling the cooling length of 
postshock gas whilst also capturing the binary orbit is a com- 
putationally demanding task and may affect the accuracy of the 
derived X-ray spectrum. Considering the cooling length of post- 
shock WR wind at periastron is ~ 6 x 10 12 cm, and that the 
finest cell size is 6 x 10 10 cm, there should be reasonably good 
sampling. However, the large temperature and density gradient 
present in the region of postshock gas in models A and B will 
result in some numerical heat conduction wh ich, as discussed 
in §|3] may affect the derived X-ray spectrum (iParkin & Pittardl 
l2010h . This situation could be improved in future models by 
adopting a higher simulation resolution in regions of large tem- 
perature and density gradients. 

5. Conclusions 

The massive WR+O-star binary system WR 22 has been investi- 
gated using 3D AMR hydrodynamic simulations which include 
radiative driving, gravity, optically-thin radiative cooling, and 
orbital motion. Two simulations were performed: one with in- 
stantaneously accelerated winds (model A), and another with ra- 
diatively driven winds (model B). We find that when the stel- 
lar winds are assumed to be instantaneously accelerated a stable 
WCR is established throughout the orbit. In this case the model 
massively over -predi cts the observed X-ray flux, consistent with 
the findings of lG09l However, in the second simulation which 
considers the acceleration of the winds, the postshock O star's 
wind transits from quasi-adiabatic at apastron to highly radiative 
at periastron and the character of the WCR changes dramatically. 
Essentially, as the stars approach periastron the ram pressure of 
the WR wind increasingly overwhelms the O star's wind, push- 
ing the WCR deeper into the O star's wind acceleration region, 
and triggering radiative cooling in its postshock wind. The sub- 
sequent growth of powerful NTSIs which massively disrupt the 
WCR is followed by a collapse of the WCR onto the O star be- 
tween ^ 0.95 - 1.05, which reduces the over-estimate of the 
observed flux by the model to a factor of ~ 6. However, dis- 
crepancies still remain: the observed flux is still overestimated 
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at all orbital phases, and the match between the models and an 
XMM-Newton spectrum at = 0.994 is poor. The latter may be 
indicating the necessity to consider in future models reductions 
in the adopted mass-loss rates, reductions in the unshocked wind 
opacity due to ionization effects, and the possibility that radiative 
braking is effective at periastron. 

Revisions to the adopted model parameters are considered 
which might improve the agreement with observations. We con- 
clude that the ratio of the mass-loss rates should be increased 
in favour of the WR star to the extent that a normal wind ram 
pressure balance does not occur at any orbital phase, potentially 
leading to a sustained collapse of the WCR onto the O star. This 
can be achieved with a ratio of WR to O star mass-loss rates 
of ^ 500, i.e. a factor of ten higher than adopted in this work. 
If radiative braking is more effective than in this work there is 
the interesting prospect of a WCR established at apastron by ra- 
diative braking transitioning to an instability-driven disruption, 
or collapse of the WCR onto the O star, at periastron. As such, 
radiative braking may play a significant role for the WCR dy- 
namics and resulting X-ray emission. 
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